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Abstract: The work is devoted to mathematical modeling of the drying process of potash fertilizers in 
fluidized bed apparatuses. The model is based on the material and thermal balance and takes into account the 
kinetics of the dehydration process of potassium chloride particles. The model is investigated using the tools of 
the Matlab mathematical computing package. The adequacy of the proposed analytical model to real processes 
and phenomena has been verified. 
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There is a significant amount of work devoted to modeling the drying processes of bulk 
materials, in particular, drying in fluidized bed apparatuses. At the same time, numerous models 
describing the structure of the fluidized bed reflect only certain aspects of the phenomenon. Due to the 
complexity of the system and with a large number of internal and external parameters, these models 
are suitable for use as a basis for reflecting empirical correlations and dependencies [1-2]. The 
mathematical model is based on the material and energy balance of drying, using simplifications and 
empirical dependencies defined, for example, in [3]. The material and energy balance of dehydration is 
determined by the following equations: 


M(1-U„)=(M-U)( 1-U;) (1) 
U =R(Ē); (2) 
qı 
qı = (quen + 2,01t,) + (m-NCyot., 5 (3) 
ma U a 
(7-0) 
q=(Q+0, 61,0) (5) 


P 


where Ux — final humidity, (mass fraction); U,, — initial humidity, (mass fraction); U— the 
amount of evaporated moisture, (kg/s), g— the amount of heat, transmitted in a layer by a coolant per 
one kg of burned fuel, (J/k); qı— the amount of heat per evaporation of one kg of moisture, (J/kg); m— 
consumption of wet material corresponding to one kilogram of evaporated moisture, (ke/ke Bularn); 
Cine — heat capacity of the material (potassium chloride-KCl), (Jkg:K); te, — fluidized bed temperature, 
(°C); t.— coolant temperature, (°C); Q- calorific value of gas, (lo1c/k2); qucn — the heat of evaporation 
of water at 0°C, (J/kg); n— coefficient of heat loss in the furnace and fluidized bed apparatus; M- the 

Published under an exclusive license by open access journals under Volume: 2 Issue: 8 in Aug-2022 


Copyright (c) 2022 Author (s). This is an open-access article distributed under the terms of Creative Commons Attribution 
License (CC BY).To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/ 


6 


eee 
J | AET International Journal of Innovative 


Analyses and Emerging Technology 
| e-ISSN: 2792-4025 | http://openaccessjournals.eu | Volume: 2 Issue: 8 


amount of wet material entering the drying, (kg/h); R- total fuel consumption, (kg/s). 


The gas temperature in the sublattice space , depending on the fuel gas consumption , is 
calculated by the empirical formula [4]: 


tr = 2920 -B / a+ 0,43. (6) 
In this case , the coefficient a is calculated by the formula 
a=L/(R-9,4) (7) 


where p — coefficient that takes into account heat loss in the furnace; a — excess air ratio; L — 
hourly secondary air consumption, ( nm’/s). 


As a result of the transformations we get 


(“2,014 (m-1)-C,,)-0,6+ (2) 
fa = a 1.2 ae 
eG 
tp 
U Q 0,6 U (8) 
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The dependence of the final humidity on the temperature of the layer is expressed by the 
empirical formula: 


Ua Ae (9) 


where the coefficients A and k depend on a variety of process conditions, in particular on the 
ambient temperature, the granulometric composition of the product, and are determined empirically 


and take the following values: A=0.02237, wt, ; k=-0,025, a 


The model based on the material-thermal balance is implemented using the means of the 
MATLAB mathematical computing package. The operation of the model was tested on an array of 
historical data from the automated process control system archive obtained on functioning equipment. 
At the same time, sampling to determine the moisture content of the product in the laboratory was 
carried out at intervals of 600 seconds during the working shift. Figure 1 shows graphs of changes in 
the humidity of potassium chloride particles at the outlet of the fluidized bed dryer. In the graph shown 
in Fig.1, the series indicated by the dotted line is a sample of the values of the final moisture content of 
the U; product in the real process (Fig. 2, a). The series indicated by a solid line represents the result of 
the analytical model at the values of its parameters indicated above, as well as the values of the 
crystallizate flow rate at the inlet and the initial humidity indicated in Fig. 2, b. If there is a selective 
distribution law of random observations (from the experiment) and the distribution law of the general 
population (determined by the model), then the adequacy of the analytical model can be assessed using 
the Pearson agreement criterion (x - criteria) [4]. 


Published under an exclusive license by open access journals under Volume: 2 Issue: 8 in Aug-2022 
Copyright (c) 2022 Author (s). This is an open-access article distributed under the terms of Creative Commons Attribution 
License (CC BY).To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/ 


7 


ee ee Se 
J | AET International Journal of Innovative 


Analyses and Emerging Technology 
| e-ISSN: 2792-4025 | http://openaccessjournals.eu | Volume: 2 Issue: 8 
To apply the y’-criterion, the entire range of variation of a random variable in a sample of 


volume N is divided into r intervals. The number of intervals when calculating the Pearson consensus 
criterion is determined by the formula 


r=14+3,322-Ig(N) (10) 
where N is the number of values in the sample, r is the number of intervals. 


The number of intervals r (for n = 45) is 6. Based on the maximum and minimum values of the 
residual humidity index in the sample under consideration and the specified number of intervals, the 
interval step of 0.02 wt is calculated. %. 


Change in humidity at the outlet 


Moisture content at the outlet of the furnace CS, 


fei te cate pose sees esi oats doe esoceceees 


Time t*10%,s 


Fig.1. Change in the humidity of potassium chloride (KCl) particles at the outlet of the fluidized 
bed apparatus 


Humidity at the entrance 


Humidity of the mold at the inlet, wt .% 


Time t*10* 5S 
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The nature of the load change 


The nature of the change in the crystallizate flow rate, y.e. 


0,5 1 1,5 2 2,5 


Time t* 107s 


Fig. 2. Results of the analytical model of the fluidized bed furnace 
Calculation of the frequencies of the results presented in Fig. (2,a), is given in Table.1 


Table.1. Results of calculation of the analytical model 


Interval Interval boundary, wt.% Values of the frequency of hits 
number in the interval 
lower upper settlement laboratory 

1 0,1 0,12 2 2 

2 0,12 0,14 8 10 

3 0,14 0,16 11 6 

4 0,16 0,18 11 9 

5 0,18 0,2 10 13 

6 0,2 0,22 3 5 


The calculated value of the y’-criterion is calculated by the formula 
“(0,-Ey 

tp? = OE = p (11) 
i=l 


where X pac- calculated value of the criterion X; O, - hit rate 


in the i-th interval of calculated values according to the model; E ; — the frequency of hitting the 
i-th interval of laboratory values; r is the number of intervals. 


The calculated value of the criterion ra ras for the data in Table.1 equals 6,5. According to the 
table of critical values of X ras » at the significance level a = 0.05 and the number of degrees of 
freedom À = r- 3 = 6 — 3 = 3 , the value of %? crete is equal to 7,8. At the same time, X° ras < X crete, i.e. 
the calculated value of the criterion turned out to be less than the critical one.Thus, the hypothesis of 
equality (agreement) of frequencies is not rejected and the model is recognized as adequate. 
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To assess the quality of the developed model, you can additionally use other indicators, for 
example: 


1. Correlation coefficient. Reflects the tightness of the linear relationship between the calculated 
values of the model and the experimental ones. 


ee Zo- YOr-y) 12) 


Don Iwp 


2. Average approximation error 


(67 -=y,) 
A 


delta = -100 (13) 


3. Standard approximation Error (RMSE) 


sko _ delta = |20» (14) 


4. Maximum approximation error 


max_ delta = max(abc(y, — y,) (15) 


Here: y, — calculated value according to the model; y> experimental value; y — estimation of 


the mathematical expectation of the vector of calculated values; y- estimation of the mathematical 


expectation of the vector of experimental values; N— the number of values in the output vector of the 
model. 


The values of the quality indicators of the analytical model (12)—(15): 


correlation coefficient Ru... ceeecececccccssssessececceceeseessseeeees 0,673 
average relative CENOM vossnscrsddassanisuceannasdesanntducepeeeeictenveinans 0,006 
RMS mean square C000. cciccstidet ee edesieeeteatedieeedeeee 0,023 
FE Vee DUT EI ETTOL s oiiire rrera encetvateecpacunstunesadalesed 0,057 


Information for assessing the accuracy of the analytical model can also be provided by the errors 
given regarding the range of values of laboratory analyses in studies and the range of corresponding 
calculations based on the analytical model. The relative reduced error is calculated by the formula 


Y = 4.100% (16) 
X, 


where Y — relative reduced error of the analytical model; A — average absolute measurement 
error; X, — the normalizing value of the scale. 


The difference between the maximum and minimum values of the laboratory analysis vector, 
equal to 0.12 wt, % was used as X, and the difference between the maximum and minimum values of 
the vector of calculated values, equal to 0.1132 wt.%. 
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The average absolute error of the analytical model is equal to: 


Yı TY i 
A= 17 
7 (17) 
The relative reduced error of the model for the range of values of laboratory analyses (16) is 
numerically equal to 14.5%. The relative reduced error of the analytical model based on the range of 
calculations of the model (16) is numerically equal to 15.3%. 
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